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ABSTRACT 

Type II topoisomerases are essential enzymes that 
regulate DNA topology through a strand-passage 
mechanism. Some type II topoisomerases relax 
supercoils, unknot and decatenate DNA to below 
thermodynamic equilibrium. Several models of this 
non-equilibrium topology simplification phenom- 
enon have been proposed. The kinetic proofreading 
(KPR) model postulates that strand passage 
requires a DNA-bound topoisomerase to collide 
twice in rapid succession with a second DNA 
segment, implying a quadratic relationship 
between DNA collision frequency and relaxation 
rate. To test this model, we used a single-molecule 
assay to measure the unlinking rate as a function of 
DNA collision frequency for Escherichia coli topo- 
isomerase IV (topo IV) that displays efficient 
non-equilibrium topology simplification activity, 
and for E. coli topoisomerase III (topo III), a type lA 
topoisomerase that unlinks and unknots DNA to 
equilibrium levels. Contrary to the predictions of 
the KPR model, topo IV and topo III unlinking rates 
were linearly related to the DNA collision frequency. 
Furthermore, topo III exhibited decatenation activity 
comparable with that of topo IV, supporting 
proposed roles for topo III in DNA segregation. 
This study enables us to rule out the KPR model 
for non-equilibrium topology simplification. More 
generally, we establish an experimental approach 
to systematically control DNA collision frequency. 



INTRODUCTION 

DNA topoisomerases are ubiquitous enzymes with essen- 
tial roles in regulating DNA topology in all domains of life 
(1-5). They are categorized into two classes, type I and II, 
based on the number of DNA strands cleaved and de- 
pendence on adenosine triphosphate (ATP) hydrolysis 
(1). Type I topoisomerases introduce a transient nick in 
one strand of duplex DNA and either pass the intact 
single-stranded segment of DNA through the nick (type 
lA) or allow controlled rotation around the phosphate 
bond in the intact single strand (type IB). Type I topo- 
isomerases relax supercoiled DNA towards thermo- 
dynamic equilibrium, as they do not use ATP. Type II 
topoisomerases (topo II), which are further subdivided 
into type IIA (eukaryotic, prokaryotic and viral) and IIB 
(archaeal, some higher plants and a few bacteria) based on 
structural, biochemical and evolutionary differences, 
introduce a transient double-strand break in one 
segment of duplex DNA (the gate or G-segment) 
through which a second segment of duplex DNA (the 
transfer or T-segment) is passed (6,7). Through this 
ATP-dependent reaction, type II topoisomerases can 
relax DNA supercoils, as well as generate and remove 
Hnks and knots in or between closed DNA molecules. 
Among the type IIA topoisomerases, DNA gyrase 
displays the unique ability to negatively supercoil DNA 
and is the only type IIA topoisomerase that uses the 
energy of ATP hydrolysis to perform work on the DNA 
(1,4,8-10). The remainder of the type IIA enzymes cata- 
lyse, on average, energetically favourable topological re- 
actions, i.e. relaxing supercoiled DNA. However, 
Rybenkov et al. (11) demonstrated that type IIA topo- 
isomerases, with the exception of gyrase, simplify DNA 
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topology to steady-state concentrations of unlinked and 
unknotted molecules that are lower than those found at 
thermodynamic equilibrium and generate topoisomer dis- 
tributions of relaxed DNA that have smaller variances 
than the equilibrium distribution. As topo II A enzymatic 
activity is coupled to the energy release of ATP hydrolysis, 
no thermodynamic principles are violated; however, the 
mechanism by which topo IIA acting at the local scale 
can sense the global topology of the DNA substrate 
remains speculative. Proposed models of non-equiUbrium 
topology simplification posit different mechanisms of 
topological discrimination, i.e. how the topoisomerase 
distinguishes DNA crossings in linked, knotted and super- 
coiled DNA from those occurring randomly in or between 
DNA molecules (8,12). Two of the models currently under 
debate, the G-segment DNA-bending model and the 
hooked-juxtaposition model, assume that local geometry 
sensing by the topo IIA enzyme influences T-segment 
capture and transfer (12-18). In the G- segment-bending 
model, the degree of G-segment bending imposed by topo 
IIA dictates the degree of non-equilibrium topology sim- 
plification (12,16,19,20). However, a recent study found 
that the degree of G-segment bending by several topo 
IIAs was similar and uncorrected with their non- 
equiUbrium unHnking activity. This finding is inconsistent 
with predictions of the G-segment-bending model and 
suggests that DNA bending alone is insufficient to 
account for non-equilibrium topology simpHfication (21). 
The hooked-juxtaposition model posits that topo IIA 
preferentially unHnks juxtaposed DNA segments that are 
sharply bent towards one another (13,14,22-24). 
Simulations suggest that the specific unHnking of hooked 
juxtapositions can account for much, but not all, of the 
observed degree of topological simpHfication (22); 
however, the model has not yet been experimentally 
tested. In contrast to these geometric selection models, 
the kinetic proofreading (KPR) model, based on KPR 
mechanisms proposed by Hopfield and Ninio (25,26), 
posits that topo IIA uses 'kinetic selection' of the 
T-segment (17-19). In this model, G-segment-bound 
topo IIA undergoes a transition to an 'activated' state 
on the first colHsion with a T-segment (Figure IB). 
Coupling this activated state transition with the energy 
of ATP binding makes it irreversible. Strand passage 
occurs only if a second T-segment coHision occurs before 
decay of the activated state. Non-equiHbrium topology 
simpHfication results because strand passage depends on 
the product of the initial, activating strand coHision fre- 
quency and the second, productive strand colHsion fre- 
quency (Figure IB). Under the assumption that these 
frequencies are comparable, the unlinking rate would 
scale with the square of the colHsion frequency. This 
non-Hnear relationship between unlinking rate and 
strand colHsion frequency results in non-equilibrium sim- 
plification, as the equiHbrium simpHfication reaction 
would depend linearly on the colHsion frequency. As a 
concrete example, consider the steady-state inter- 
conversion of linked and unlinked circular DNA mol- 
ecules. At equiHbrium, the concentrations of linked and 
unHnked DNA would be proportional to the ratio of 
strand collision frequencies of unlinked and Hnked 



circles. Under the KPR model, the concentrations of the 
two species would be proportional to the square of the 
DNA colHsion frequencies for the unHnked and linked 
molecules. The quadratic dependence on the coHision fre- 
quency arises because, in order for strand passage to take 
place, two coHisions must occur in rapid succession, the 
probabiHty of which scales as the square of the coHision 
rate. Thus, the central tenet of the KPR model is that the 
strand transfer rate scales quadratically with the DNA 
strand-coHision frequency. 

Simulations of the KPR model suggest that the kinetic 
amplification may not be sufficient to achieve experimen- 
taHy observed levels of non-equiHbrium topology simpH- 
fication for DNA knots unless combined with aspects of 
the DNA G-segment-bending model or some other com- 
plementary non-equilibrium simplification mechanism 
(16,27). Furthermore, the KPR model seems difficult to 
rationalize with accepted models of strand passage in 
which ATP binding triggers N-gate closure that traps 
the T-segment and favours its passage through the 
cleaved G-segment (7,28-31). This model, based on a com- 
bination of structural and biochemical evidence, seems to 
be at odds with the central tenet of the KPR model that 
postulates T-segment release and recapture after the 
binding of at least one ATP (8,10). Despite this apparent 
contradiction, there is sufficient uncertainty in the kinetics 
of binding two ATP molecules and the subsequent kinetics 
of N-clamp closure that the KPR model, which is a kinetic 
model, cannot be entirely ruled out based on existing data. 
Therefore, the hypothesis that topo IIA requires a double 
coHision between the G-segment-bound topoisomerase 
and a T-segment for strand passage remains a potential, 
but unverified, mechanism that would result in non- 
equiHbrium topology simplification. Although ensemble 
approaches allow geometric features of DNA- 
topoisomerases complexes to be indirectly measured and, 
therefore, used to test the validity of models based on 
DNA bending, testing the KPR model experimentaHy is 
challenging, as it requires a means of assessing and 
varying the coHision frequency between DNA strands, 
referred to hereafter as colHsion frequency. 

Here, we tested the KPR model by investigating the 
relationship between DNA colHsion frequency and DNA 
unlinking rate using a novel strategy based on the 
micromanipulation of a single-DNA crossing. We 
measured the unlinking rates of Escherichia coli topo IV, 
a type IIA topoisomerase that displays efficient non- 
equiHbrium topology simplification, using a single- 
molecule DNA unlinking assay in which the relative 
colHsion frequency could be controHed (Figure lA). 
Combining these measurements with colHsion frequencies 
obtained from Monte Carlo (MC) simulations permitted a 
sensitive measure of the relationship between colHsion fre- 
quency and unlinking rate. Because this comparison reHes 
on the assumption that the colHsion frequency obtained 
from MC simulations is proportional to the actual coHi- 
sion frequency, we measured the unHnking rate as a 
function of computed coHision frequency for E. coli 
topo III, a type lA topoisomerase that unHnks and 
unknots DNA to equilibrium levels (11). Topo III is an 
efficient decatenase that is thought to play important roles 
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Figure 1. Topoisomerase unlinking assay and KPR model. (A) Cartoon of experiment (not to scale). A l-jim magnetic bead (beige) is tethered to the 
surface of the flow chamber by two 5kb torsionally unconstrained dsDNA molecules (blue and gold). A Hnk between the DNA molecules is 
generated by rotating the bead (black arrow) using a magnet assembly (not shown), which decreases the DNA extension by AZ (red arrow). Strand 
passage mediated by a topoisomerase (shown bound to DNA), unhnks the DNA (blue arrow), leading to an increase in DNA extension. After a brief 
re-equihbration time, the process is automatically repeated. Extension of the double tether as a function of time shows the decrease in extension 
associated with introducing the hnk, and the subsequent increase in extension when the hnk is resolved by the strand-passage reaction of the 
topoisomerase. The time between these events is defined as the waiting time (r^ait)- DNA extension was measured by video tracking the bead in real 
time at 100 Hz with ~2-nm spatial resolution. (B) DNA strand-passage mechanisms. DNA unhnking can occur via a single- or double-collision 
mechanism. Topoisomerase binds the G-segment (blue dot 1) and collides with a potential T-segment with rate constant (blue dot 2). In the 
single-colhsion mechanism (blue shaded arrow), strand passage occurs from this state with rate constant /Cul- In the double-colhsion mechanism, i.e. 
the KPR model (pink shaded arrow), the initial T-segment colhsion activates the topoisomerase (yellow dot 1*). The activated topoisomerase will 
effectuate strand passage if it colhdes with a second potential T-segment (yellow dot 2*) before it decays back to the inactive state (blue dot 1) with 
rate constant /cd- The single-colhsion strand-passage rate (Rul) is proportional to the colhsion frequency k^, whereas the double-colhsion Rul, 
posited by the KPR model, is proportional to the square of the collision frequency. Here, /con and /coff are on and off rate constants of topoisomerase 
binding to G-segment DNA. k^. and /c/ are release rate constants of a potential T-segment by the topoisomerase. 



in chromosome segregation and the resolution of hnked 
DNA intermediates without using ATP (32-35). Topo III 
uses a strand-passage mechanism common to all type lA 
topos (1,4,36,37), but it can uniquely pass duplex DNA 
through a single-stranded nick or gap (1), thereby 
permitting reversible linking and knotting reactions with 
gapped DNA substrates. In the context of the KPR 
model, topo III is expected to display a Hnear relationship 
between unlinking rate and colhsion frequency because it 
uses an ATP-independent mechanism based on a single 
colhsion of two DNA segments (Figure IB). 

Contrary to the predictions of the KPR model, we 
found that the unlinking rates of both topo IV and topo 
III were hnearly related to the collision frequency. 
Moreover, the unlinking rates of the two enzymes 
measured with comparable double tether geometries 
were hnearly related over the range of DNA colhsion 



frequencies investigated. We conclude from these results 
that non-equilibrium topology simplification activity does 
not arise in whole or in part from a KPR mechanism, thus 
ruhng out the KPR model. However, our results reveal 
that topo III, unlike E. coli topoisomerase lA (topo I), 
is an efficient DNA decatenase, on par with topo IV, 
lending support to its putative roles in chromosome seg- 
regation and unlinking in vivo (33). 

MATERIALS AND METHODS 

Topoisomerase IV unlinking substrate 

Linear 5-kb dsDNA labehed with a single biotin or 
digoxigenin at each 5^-end was prepared by polymerase 
chain reaction (PCR) using pET28b (EMD Science) as a 
template with a biotin-labelled forward primer and 
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digoxigenin-labelled reverse primer (Eurofins MWG 
Operon). The PGR product was purified with a 
QIAquick PGR purification kit (Qiagen). 

Topoisomerase I and III unlinking substrate 

Linear 5-kb dsDNA with a 37-nt gap (2258-2294) with a 
single biotin or digoxigenin at each 5^-end was prepared 
similarly to the topo IV DNA substrate except that the 
PGR template was pKZl. pKZl was made by Hgating a 
44-bp dsDNA segment containing two BbvGI sites 
separated by 37 nt into pET28b at the BamHI site. The 
5-kb PGR product was digested with Nt.BbvGI (New 
England Biolabs) to create two nicks at positions 2258 
and 2294, near the middle of the DNA template. After 
digestion, the 37-nt oHgomer was removed by heating 
the DNA to 80°G for 20min and cooling down gradually 
to 4°G in the presence of excess complementary 37-nt 
oligomer. The 4-kb multi-nicked DNA molecules were 
prepared by PGR of pFG94 (38) from position 5422- 
9520 with a biotin-labelled forward primer and 
digoxigenin-labelled reverse primer (Eurofins MWG 
Operon) and were nicked with four nicking enzymes: 
Nb.BsmI, Nt.AIwI, Nb.BtsI and Nt.BtsNBI 
(New England Biolabs), generating 14 nicks located over 
the middle third of the DNA. 

Protein preparations 
E. coli topoisomerase I 

E. coli topoisomerase I was purchased from New England 
Biolabs. 

E. coli topoisomerase III 

The pET-15b^''^^^ plasmid encoding K12 coli topoisom- 
erase III was provided by J. Keck (University of 
Wisconsin Madison). A tobacco etch virus (TEV) 
protease cleavage site with an additional two glycine 
residue spacer was inserted upstream of the topB gene 
using the QuikGhange II XL mutagenesis kit 
(Stratagene), and the sequence of the resulting plasmid 
(pET-15b'^^^+^^^) was verified. Rosetta (DE3) pLysS 
Gompetent Gells (EMD4Bioscience) transformed with 
pET-15b'^^^+^^^ were grown at 37°G in Luria-Bertani 
broth containing 50|ig/ml of carbenicillin and 34|ig/ml 
of chloramphenicol. The expressed protein was induced 
with 1 mM isopropyl P-D-thiogalactopyranoside. After in- 
duction at 37°G for 20 h, cells were harvested by centrifu- 
gation. The collected cells were suspended in buffer A 
[50 mM Tris-HGl, pH 7.5, 300 mM NaGl and 10% (v/v) 
glycerol] and emulsified using a high pressure cell hom- 
ogenizer (Avestin). The cell lysate was then centrifuged at 
40krpm at 4°G for 1 h in a Ti45 rotor (Beckman). The 
supernatant was loaded onto a HisTrap-FF 5 ml 
Ni-Nitrilotriacetic acid (NT A) column (GE). His- tagged 
protein was eluted with buffer B [50 mM Tris-HGl, pH 
7.5, 0.5 M imidazole, 0.3 M NaGl, 10 mM 
P-mercaptoethanol and 10% (v/v) glycerol]. The eluted 
proteins were concentrated and buffer-exchanged with 
buffer A using an Amicon Ultra- 15 centrifugal filter unit 
with Ultracel-30 membrane (Millipore). The his tag was 
removed from the protein by digestion with 0.5 U/|ig of 



Pro-TEV (Promega) for 6 h. The protein digestion mixture 
was run twice through an Ni-NTA column (Qiagen) to 
remove Pro-TEV and uncleaved protein. The size and 
purity of cleaved topo III were evaluated by sodium 
dodecyl sulphate-polyacrylamide gel electrophoresis 
(Invitrogen). 

E. coli topoisomerase IV 

ParG and ParE subunits of topo IV were purified to a final 
concentration of 1.5 |iM as previously described (21,39). 
Equimolar amounts of ParG and ParE were mixed to 
make 0.75 |iM of topo IV heterotetramer. 

Experimental procedure 

To prepare double-DNA tethers for DNA unlinking ex- 
periments, 1 |il of 0.15nM singly biotin- and digoxigenin- 
labelled 5-kb dsDNA was incubated with 1 |il of magnetic 
beads [1% (w/v), MyOne, Invitrogen] in 200 |il of wash 
buffer [WB: Ix phosphate-buffered saline, 0.3% (w/v) 
bovine serum albumin (BSA) and 0.04% Tween-20) at 
room temperature (22° G) overnight at a 3:1 DNA to 
bead ratio. The sample cell was incubated with 30 |il of 
10|ig/ml of anti-digoxigenin in Ix phosphate-buffered 
saHne for 1 h at room temperature followed by a wash 
with 200 |il of WB. The DNA-bead mixture was incubated 
in the anti-digoxigenin-coated sample cell for 30min 
before washing with 400 |il of WB. The magnetic 
tweezers instrumentation has been previously described 
(40-42). During the experiment, double-DNA tethered 
beads were identified by rotating the magnet assembly. 
In contrast to a bead tethered by a single DNA, the 
DNA extension of a double-DNA tethered bead decreases 
symmetrically for negative and positive turns (Figure 2A) 
independent of the tension on the DNA (43). Topo IV 
activity buffer [25 mM Tris-HGl, pH 7.5, 100 mM potas- 
sium glutamate, lOmM magnesium chloride, 1 mM 
dithiothreitol, 0.3% (w/v) BSA and 0.04% (v/v) 
Tween-20] containing 2 nM topo IV and 1 mM ATP was 
introduced once 20 double-tethered beads were found, and 
their tethering geometries (spacing between DNA mol- 
ecules and verification of parallel tether geometry — see 
later in the text) were characterized. Topo IV DNA un- 
linking rates were measured for bead rotations of —0.6 to 
— 1 turn (Figure 3 A), which imposes left-handed DNA 
crossings as found in positive writhe. Topo III unHnking 
measurements were similar except that 5-kb DNA mol- 
ecules containing a 37-nt gap near their centres were 
used as substrates, as topo III requires a single-stranded 
region for strand passage. Unlinking measurements were 
performed with 5nM topo III in topo III activity buffer 
[40 mM HEPES, pH 8.0 adjusted with KOH, 1 mM mag- 
nesium acetate, pH 7, 0.3% (w/v) BSA and 0.04% (v/v) 
Tween-20]. Data traces were analysed using a custom 
step-finding algorithm to identify the waiting time, T^ait 
(the time interval between rotating the bead to generate 
the DNA Hnk and it being unHnked by topoisomerase). 
The Twait values (100-200 events for each rotation) were 
binned, and the histogram was fitted with a single expo- 
nential to obtain the strand-passage rate. Topo I (New 
England Biolabs) unHnking measurements were 
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Figure 2. Characterization of DNA tether geometry and 
rotation-dependent coUision frequency obtained from MC simulations. 
(A) In situ measurement of DNA tether geometry (not to scale). The 
double tether consists of two 5-kb DNA molecules attached between 
the surface and a streptavidin-coated magnetic bead (grey). To obtain 
the tether geometry, the DNA extension was measured as the bead was 
rotated from —0.5 to 0.5 rotations in increments of 0.1 turns (blue dots 
and dashed Hne). The DNA extension over the range —0.5-0.5 rota- 
tions was fit with Equation (1) (black sohd line) to obtain the spacing 
between the DNA molecules, 2e, and the DNA extension at zero 
rotation, Zq. (B) Example colhsion frequency distributions for different 
imposed turns, n, from MC simulations of one double tether. 
Distributions of colhsion frequencies (number per 10^ simulation 
steps) are fit with Gaussians (black sohd hues) to obtain the mean 
colhsion frequency, <fsim>, for each turn. (C) Mean colhsion fre- 
quency, <fsim>, as a function of imposed turns, n. Values are 
colour-coded based on the corresponding number of imposed turns 
shown in B. 

performed using the same substrate used for topo III ex- 
periments with ~ 130 or 325 nM of topo I in topo I activity 
buffer [50 mM potassium acetate, pH 7.9, lOmM magne- 
sium acetate, 1 mM dithiothreitol, 0.3% (w/v) BSA and 
0.04% (v/v) Tween-20]. 



DNA crossing geometry calculation 

The geometry of the two DNA strands tethering the bead 
to the surface was determined by measuring the extension, 
Z(n), as a function of the number of turns (n). The region 
from —0.5 to 0.5 turns was fit with Equation (1) to obtain 
the distance between two DNA molecules (2e) and the 
DNA extension at ^ = 0 (Zq) (Figure 2A): 

Z(n) = ^Zl - Ae^sm^ilnn) - Rb^^^l (1) 

where is the bead radius (0.5 |im) (44). Tethered beads 
were selected for unlinking measurements if the DNA ex- 
tension at zero rotation (Zq) was within 10% of the pre- 
dicted extension at the applied force. This selection 
criterion eliminated tether geometries that were not 
parallel, as trapezoidal tether geometry would result in a 
decreased DNA extension at zero rotation, and eliminated 
beads tethered by more than two DNA molecules. 

DNA collision frequency calculation from Monte Carlo 
simulations 

The experimentally determined tether geometry (2e, n, 
DNA length and applied force, F) was used to calculate 
the number of coUisions between two DNA segments 
using a custom-written MC simulation routine for 
double-tethered beads, as previously described (44,45). 
Briefly, the simulation program models the two DNA 
molecules as discrete worm-Hke chains consisting of 10- 
nm rigid segments of fixed length. For each braid 
geometry, 10 simulation steps were calculated and were 
sampled once every 10^ steps to obtain uncorrected con- 
figurations. At each sample step, the colhsion frequency, 
fsim, was calculated from the total number of colHsions 
over the previous 10^ MC steps. ColHsions were con- 
sidered to have occurred when two segments on the two 
DNA molecules were closer than a cut-off distance of 
12 nm. Thus, for 10^ steps, 10^ colhsion frequencies were 
generated for each tether geometry. The fsim values for 
each condition were binned, plotted as a histogram and 
fitted with a Gaussian distribution to obtain the mean fsim 
as shown in Figure 2B. 

RESULTS 

Single-DNA crossing assay measures unlinking rate as a 
function of DNA strand collision frequency 

To test the KPR model of non-equihbrium topology sim- 
pHfication, we developed a single-molecule DNA unlink- 
ing assay (Figure 2A) that permitted measurement of the 
unhnking rate as a function of strand colhsion frequency. 
The KPR model predicts a quadratic relationship between 
unhnking rate and colhsion frequency. The unlinking 
assay is a variation on similar unlinking assays (44-46) 
in which a magnetic bead is tethered by two rotationally 
unconstrained 5-kb DNA molecules. Rotating the bead 
wraps the two DNA strands around one another, which 
reduces the height of the bead. The action of a single 
topoisomerase on this substrate resulted in a sudden 
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Figure 3. DNA extension traces of unlinking and mean rates of unlinking as a function of bead rotation. (A) DNA extension traces show unlinking 
activity of topo IV at three different imposed bead rotations (0.6, 0.8 and 1) for one double tether. The distributions of waiting times (Twait) (right 
panels) are well fit by single exponentials (black soHd fines) to obtain tfie unfinking rates, k(n). (B) Tfie unfinking rates are plotted as a function of 
bead rotation, n. Error bars correspond to tfie SD of tfie fit parameter. 



increase in the height, indicating the unhnking of the 
DNA strands. After a brief delay, the bead was automat- 
ically rotated to create a new substrate for the enzyme 
(Figure 1). The experiment consists of imposing a series 
of partial rotations (— 1</7<— 0.6, corresponding to 
positive writhe) on the bead and measuring the rate at 
which the corresponding crossings are unHnked by topo 
IV. The frequency of DNA strand colHsions is expected to 
increase in proportion to the degree of Unking, i.e. the 
magnitude of the bead rotation, which has been verified 
by MC simulations (Figure 2B and C). Therefore, by 
measuring the unhnking rate as a function of imposed 
bead rotation and estimating the strand colHsion fre- 
quency from MC simulations of the experimental config- 
uration, the relationship between strand colHsion 
frequency and unlinking rate can be determined. 

In practice, once a double-tethered magnetic bead 
was identified, the DNA crossing geometry was 
characterized to obtain the spacing between the DNA 
tethers, 2e (Figure 2A). Based on the tether geometry 
and the applied force, the inter-strand colHsion frequency 
as a function of bead rotation was obtained from MC 
simulations (Figure 2B and C). The unhnking rate was 
measured as a function of bead rotation from —0.6 to 
— 1 turn in the presence of 2nM topo IV (above the 
of ~0.2nM measured under similar conditions) (47) and 



saturating (1 mM) ATP. Under these conditions, 
T-segment capture is expected to be rate limiting for 
strand passage; thus, the unlinking rate will be determined 
by the DNA collision frequency. The observed increase in 
unlinking rate with rotation (Figure 3B) is consistent with 
this assumption. The time interval between rotating the 
bead and the unlinking reaction was obtained from the 
extension traces (Figure 1), and the distributions of un- 
linking times were fit with single-exponential functions 
to obtain the mean unlinking rate for each rotation 
(Figure 3). To compare unhnking rates as a function of 
bead rotation measured with different double-tether 
geometries and to account for other possible differences 
among nominally identical measurement conditions, un- 
linking rates for each tether were normalized to the un- 
linking rate at ^ = 1 turn (Figure 4A). 

The unlinking rate of topo IV is linearly related to the 
collision frequency 

The normalized unhnking rates, kr-iv, increased as a 
function of bead rotation at both 1 and 2 pN of applied 
force on the bead (Figure 4A). To investigate the relation- 
ship between the unhnking rates and colHsion frequencies, 
we plotted kr-ivin) as a function of the normalized 
collision frequencies at the same bead rotation, fsimin) 
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Figure 4. Normalized unlinking rates of topo IV as a function of bead 
rotation and normalized collision frequency. (A) Normalized unlinking 
rates of topo IV measured with 12 different double tethers, kr-iv, are 
plotted as a function of bead rotation (n) for apphed forces of 1 pN 
(triangles) and 2 pN (circles). Error bars are calculated through propa- 
gation of uncertainties (SD) in the initial fit parameters. (B)_kT-iv 
as a function of the normalized simulated colHsion frequency, fit 
with Hnear (soHd fine), quadratic (dashed fine) and power (short-long 
dashed fine) functions. The power fit, kr-iv = Afsim returned a power 
close to 1 (^ = 0.90 ± 0.01 and P = 1.10 ± 0.03, xl = 9.5). The data 
are poorly fit by a quadratic relation (kr-iv = ^f^im) evidenced by 
the 2-fold increase in reduced chi-square, xl = 19.2. Fitting to a fine 
returned a slope of 0.87 ± 0.01 and a Hnear correlation coefficient of 
0.75 with xl = 9.5. 



(Figure 4B). The fsim values were_ calculated in a similar 
manner to that used to obtain kr-iv' the collision fre- 
quency at each rotation was normahzed by the colli- 
sion frequency at /7 = 1 turn for each double-tether 
geometry. The normalized unHnking rates seemed to be 
linearly related to the normalized collision frequencies 
(Figure 4B). Fitting^ kr-iv versus fsim with a power 
function, kj-iv = ^f^uiv returned a power close to 1: 
P= 1.10 ±0.03 (±SD) with a reduced chi-square, 
xl = 9.5. Fitting with P = 2 resulted in a 2-fold increase 
in reduced chi-square, = 19.2. Fitting to a line {P = 1) 
returned a slope of 0.87 ± 0.01 (±SD) and a linear cor- 
relation coefficient of 0.75. Overall, the fits suggest a linear 
relationship between the normalized rates and colHsion 
frequencies in contrast to the quadratic relationship pre- 
dicted by the KPR model. 



Topo III but not topo I efficiently decatenates linked 
DNA molecules containing a single-stranded gapped 
region 

Although we found a linear relationship between the 
normalized unlinking rates of topo IV and the predicted 
colHsion frequencies, it is possible that the colHsion 
frequencies obtained from MC simulations were not 
good estimates of the relative colHsion frequencies. To 
verify the MC simulation results, we measured the unlink- 
ing rates of topo III as a function of bead rotation and 
compared the normalized rates with the normalized 
computed coHision frequencies. 

The unlinking rates of topo III (10 nM) as a function of 
bead rotation were measured using the same protocol and 
conditions used for topo IV except that the 5-kb DNA 
tethers contained a 37-nt single-stranded gapped region 
near the middle of each DNA molecule, as intact 5-kb 
DNA molecules and multi-nicked 4-kb DNA molecules 
were not decatenated by topo III. These observations con- 
firmed that a single-stranded region of DNA is essential 
for topo III decatenation in addition to relaxation, similar 
to topo I (35,36). However, we found that E. coli topo lA 
was incapable of decatenation. Despite rapidly relaxing 
negative supercoils, topo lA did not unHnk crossings of 
a double-DNA tether at a ~30-fold higher concentration 
(325 nM) than that of topo III, for which efficient unlink- 
ing of the same tether was observed (Supplementary 
Figures S2 and S3). The average unHnking rate at one 
turn for topo III (5.9 ± 0.7 s~^ mean ± standard error) 
was comparable with that for topo IV (5.1 ± 0.5 s~^), 
indicating that topo III is an efficient decatenase on par 
with topo IV. The normalized_ rates for topo III (kr-iii) 
were plotted as a function of fim, calculated in the same 
manner as for topo IV (Figure 5A) and fitted with a Hnear 
function. As expected, the topo III unlinking rate was 
linearly related to, and strongly correlated with, the coHi- 
sion frequency. The fit returned a slope = 1.00 ± 0.01 
(±SD), = 6.8, and a linear correlation coeffi- 
cient = 0.76. 

Topo III and topo IV display similar rotation-dependent 
unlinking rates 

Comparing the unlinking rates of topo III and topo IV as 
a function of bead rotation, i.e. DNA coHision frequency 
for similar bead geometries, provides an independent test 
of the KPR model that does not depend on the results of 
the MC simulations. As the unHnking rate of topo III is 
reasonably expected to be linearly related to the strand 
colHsion frequency, the normalized topo III unlinking 
rate reflects the relative strand coHision frequency. 
Therefore, plotting the normalized relaxation rates of 
topo IV as a function of the normalized relaxation rates 
of topo III at the same bead rotation for comparable 
tether geometries (difference in 2^<0.2|im) and applied 
force provides a simulation-independent comparison 
between unlinking rate and strand coHision frequency 
(Figure 5B). Fitting kr-iv versus kr-iii to a line returned 
a slope of 0.97 ± 0.02 with xl = 5.8 and a linear correl- 
ation coefficient of 0.78. Fitting the data with a power 
function, kr-iv = A(kT-iiif, returned A = 0.99 ± 0.04 
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Figure 5. Normalized unlinking rates of topo III versus normalized 
collision frequency and comj^arison with topo IV. (A) Normalized un- 
linking rates of topo III, kj_ui, are plotted as a function of the 
normalized collision frequency from MC simulations, fsim, and fit to 
a line, with slope = 1.00 ± 0.01, Xv = 6.8, and a linear correlation co- 
efficient = 0.76. Error bars are calculated through propagation of 
uncertainties (SD) in the initial fit parameters. (B) Normalized rates 
of topo IV are directly compared with the normalized rates of topo 
III in a scatter plot, in which each point corresponds to the two 
normalized unlinking rates (topo III -x, and topo IV -y) measured 
with double tethers of comparable geometry at the same imposed 
rotation. Fitting the scatter plot with a Hnear function gave a slope 
S = 0.97 ± 0.02 with linear correlation coefficient 0.78 and xl = 5.8. 
Fitting to power function (y = Ax^) returned A = 0.99 ± 0.04 and 
P = 1.02 ± 0.05 with xl = 5.2. Fitting to a quadratic function 
(y = Ax^) resulted in a significantly increased reduced chi-squared 
value, Xv = 



and P = 1.02 ± 0.05 with xl = 5.2. Fixing P = 2, i.e. a 
quadratic fit, increased xl to 17.9. This direct comparison 
confirmed that the strand-passage reaction of topo IV and 
topo III depends Hnearly on the DNA colHsion frequency. 



CONCLUSION AND DISCUSSION 

In this study, we showed that the unHnking rate of topo IV 
increases linearly rather than quadratically with increasing 
colHsion frequencies between two DNA segments. The 
linear relationship between the topo III unlinking rate 
and colHsion frequency confirmed the vaHdity of the 
combined experimental and computational approach 
while providing a simulation-independent test of the 
KPR model. OveraH, our study indicates that a single col- 
lision with T-segment DNA is likely sufficient for strand 
passage by topo IV, suggesting that non-equiHbrium 



topology simpHfication cannot be explained within 
the framework of the KPR model (17,18). This conclusion 
is further supported by the fact that Drosophila topo II 
(48) showed a similar Hnear relationship when compared 
with topo III under similar geometric conditions 
(Supplementary Figure SI). Although theoretical and 
computational studies have shown that the KPR model 
could not fuHy explain the measured non-equilibrium 
topology simplification (12), our work experimentaHy 
demonstrates that topo IV does not require a double col- 
Hsion with a T-segment DNA for strand passage. This 
result Hkely applies to other topo IIA enzymes, as they 
share a common two-gate mechanism for a strand 
passage through weH-conserved core domains (4,49,50), 
and topo IV exhibits the most dramatic non-equilibrium 
topology simpHfication activity among the type IIA 
topoisomerases for which this activity has been 
characterized (11). 

Since the remarkable discovery by Rybenkov et al. (11) 
that type IIA topoisomerases can simplify DNA topology 
to below equilibrium levels, several mechanisms of 
non-equiHbrium topology simplification have been 
proposed, and their merits have been debated, largely 
through simulations. Recently, experimental tests of 
some non-equiHbrium models have been performed. For 
example, direct measurements of the bend angles imposed 
on G-segment DNA by three type IIA topoisomerases 
suggest that the bend-angle model can account for some 
but not all of the non-equilibrium simplification activity 
(21). Stuchinskaya et al. (27) made a number of measure- 
ments of non-equiHbrium topology simpHfication activity, 
some of which ruled out a 'translocation' model originally 
proposed by Rybenkov et al. to explain the non- 
equiHbrium activity. In this model, type II topos use the 
energy of ATP hydrolysis to translocate along the DNA, 
thereby localizing Hnks to smaHer loops that are more 
readily removed. Contrary to the premise of this model, 
the presence of high affinity DNA-binding proteins acting 
as translocation 'roadblocks' did not alter the degree of 
non-equiHbrium topology simplification. A 'three- 
segment' binding model that postulates the binding of 
two separate potential T-segments, one of which is 
passed through the G-segment, provides a mechanism of 
non-equiHbrium simplification, but it is also associated 
with a shift of the centre of the steady-state topoisomer 
distribution that has not been observed in multiple inde- 
pendent studies (21,27). The remaining models that have 
not been directly tested are the KPR model and the 
hooked-juxtaposition model. The hooked-juxtaposition 
model shares elements of the bend-angle model, but 
rather than the enzyme imposing a bend on the 
G-segment, the model postulates that the enzyme prefer- 
entially unHnks T- and G-segments that are bent towards 
each other. We focused our attention on the KPR model, 
in which the strand-passage rate scales non-linearly, likely 
quadraticaHy, with the strand coHision frequency. 

The experimental challenge was to devise a method of 
reproducibly controlling the strand coHision frequency 
and measuring the unHnking rate. To address this chal- 
lenge, we modified a single-molecule DNA 'braiding' 
assay, in which two strands of DNA attached to a 
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magnetic bead can be wrapped around one another by 
rotating the bead in a magnetic tweezers instrument. We 
reasoned that rotating the bead by less than one full turn 
would impose on average a partial link between the two 
DNA strands, and that the relative frequency of strand 
collisions that a topoisomerase could act on would 
increase with increasing bead rotation. MC simulations 
of the DNA configurations based on the geometry of the 
double tether measured in situ confirm that the relative 
probabiHty of strand collisions, which is proportional to 
the frequency of strand coUisions under reasonable equi- 
librium assumptions, increases Hnearly with increasing 
bead rotation (Figure 3). Unlinking measurements with 
topo III confirm the increase in colHsion frequency with 
increasing bead rotation and provide a test of the assump- 
tions impHcit in comparing the results of MC simulations 
that provide the probabiHty of strand colHsion with ex- 
perimental measurement of unlinking rate. As the conclu- 
sions of the current study rely heavily on this reasonable 
but untested hypothesis, the topo III results were an im- 
portant verification of the approach and serve as a valid- 
ation of the method, which may be generally useful to 
probe protein-mediated DNA synapse formation. In 
addition, the unlinking measurements of topo III 
provide explicit experimental evidence, indicating that 
topo III is an efficient decatenase, comparable with topo 
IV, provided the Hnked DNA includes a single-stranded 
gapped region, lending support to its putative roles in 
chromosome segregation and unlinking in vivo (33). In 
comparison with topo III, we found that topo I failed to 
unHnk crossed DNA molecules when tested with the same 
substrate, confirming that topo I and topo III likely play 
different roles in in vivo (Supplementary Figure S3). 

The combination of the partial linking assay and MC 
simulations of DNA conformations permitted a sensitive 
test of the non-linear relationship between strand colHsion 
frequency and unlinking rate that is the central feature of 
the KPR model. Ruling out the kinetic selection embodied 
in this model raises the question as to how topo IIA 
enzymes achieve non-equilibrium topology simpHfication. 
As has been pointed out by Vologodskn and others, local 
DNA conformations and the corresponding probabiHty of 
topological conversions by a topoisomerase Hkely result in 
different degrees of non-equilibrium topology simplifica- 
tion by different topo IIA enzymes (11,13-17,19). As 
proposed models have been experimentaHy tested and 
shown to be either inconsistent or insufficient to fully 
account for the degree of topological simpHfication, we 
posit that either local DNA bending coupled with an as 
yet undetermined mechanism or a selection based on an 
existing local conformation such as 'hookedness' hold the 
greatest promise to explain this phenomenon (13,14). MC 
simulations of the DNA tethers used in this study suggest 
that local DNA bending and hookedness did not vary 
significantly under our experimental conditions; thus, we 
were unable to test the effects of local geometry with the 
current measurements. We anticipate that similar 
approaches may be applicable to evaluate the effect of 
local bending and hookedness on DNA unlinking 
activity, and that these approaches will help elucidate 
the mechanistic basis for non-equilibrium topology 



simplification and, more generally, the process of strand 
capture and passage that is central to the activity of type II 
topoisomerases. 
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Supplementary Data are available at NAR Online: 
Supplementary Figures 1-3 and Supplementary 
References [48,51]. 
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